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i 1 Abstract 

Flat directions in the Minimal Supersymmetric Standard Model (MSSM) can become unsta- 
ble due to radiative corrections, and the global minimum of the (zero temperature) potential 
can lie at large values of the squark and/or slepton fields. Here we show that, in inflation- 
ary models of early cosmology, the universe is very likely to be in the domain of attraction 
of this global minimum at the end of inflation. While the minimum at the origin of field 
space is global at sufficiently high temperatures, depending on details of the model, the 
universe may be trapped in the non-zero minimum until it becomes the global minimum at 
low temperatures. Parameter values leading to this scenario are therefore ruled out. 



In an earlier paper we pointed out that certain mass patterns for the superpartners 
result in a radiatively-corrected scalar potential whose global minimum is at large values 
of the squark and/or slepton fields. Because these minima are in general electric charge 
and color breaking, we concluded that these mass patterns were therefore disallowed, or 
at least required new physics beyond the standard model. However, Riotto and Roulet 
argued that the desired minimum (with zero squark and slepton VEVs and the usual nonzero 
Higgs VEVs), while not global at zero temperature, is generically the global minimum at 
high temperature. Furthermore, if, for some reason, the universe begins in the domain of 
attraction of this minimum, then it will get stuck there, since the probability to tunnel to 
the global minimum (after the temperature drops sufficiently) is much too small. A similar 
argument was made in Q. 

The assumption made (either explicitly || or implicitly 0) is that the universe will begin 
in the global minimum of the temperature dependent scalar potential. This, however, is 
extremely unlikely to be the case along flat directions in field space after a period of inflation 
Q. During inflation, fluctuations in scalar fields which parameterize flat directions, i.e. those 
with supersymmetry breaking squared masses m 2 <C H 2 increase with time according to || 

(<P 2 ) = ^jflV , (1) 

where H is the Hubble constant during the inflationary epoch. Subsequently, the long 
wavelength modes of these fluctuations behave as a classical background field with amplitude 
ip ~ (ip 2 ) 1 ^ 2 . More precisely, the value of (po we compute in this way is the RMS width 
of a probability distribution (which is approximately Gaussian) for the value of if in a 
horizon volume. Solving the horizon and flatness problems requires that the duration r of 
the inflationary era satisfies Ht > 60. In very general models of inflation which involve a 
single mass scale mj, we have Hr ~ (Mp/mi) 2 , where Mp is the Planck mass. The value of 
raj determines the size of the observed fluctuations in the microwave background, and their 
measured value fixes (mi/Mp) 2 ~ few x 1CT 8 @, which then implies Hr > 10 8 [[/]]. The 
case of interest for us will be that of a squark or slepton field (p (which is not the inflaton 
field) with a slowly varying potential V(ip) corresponding to an almost flat direction. We 
define the field-dependent squared mass m 2 ((p>) = (1 / 2)d 2 V \ip) / dip 2 . The field stops evolving 
according to eq.(|l]) when m 2 (ip)t ~ 3H/2. Then we find |J 

For H 2 ~ mj/Mp ~ 10 _15 Mj,, and for (nearly) flat directions characterized by a mass scale 
m(ip) ~ 10~ 16 Mp, the first term in eq. (0) is smaller than the second, and the growth of 
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fluctuations is therefore limited by the duration r of the inflationary era. This generally 
results in a value of (po which is much larger than H, tpo ~ 10~ 4 Mp. 

Thus, if ifo is as large or larger than the location of the global minimum of V(ip), we 
are overwhelmingly likely to find ourselves, after inflation, in the domain of attraction of 
this minimum. This fact (that scalar fields along flat directions find themselves far from the 
origin in field space after inflation) has long been recognized and is the basis of the Affleck- 
Dine mechanism of baryogenesis || ||, which requires large initial values of squark and 
slepton fields. This scenario is particularly attractive because one can in general find flat 
directions in field space in supersymmetric grand unified theories. That is, these theories 
naturally contain directions in which the scalar potential is absolutely flat V((p) = up to 
supersymmetry breaking effects which induce masses of order the supersymmetry breaking 
scale, ip is some combination of squark and slepton and Higgs fields in which the F- and D- 
terms in the scalar potential vanish. If the expectation value of <p is non-zero (and large) 
initially (e.g. due to inflation as described above), the subsequent evolution of certain flat 
direction in a supersymmetric GUT can be shown to give rise to a baryon asymmetry. 

It has recently been pointed out however, that the simple picture of driving scalar fields 
to large vacuum values along flat directions during inflation is dramatically altered in the 
context of supergravity flTU|. During inflation, the Universe is dominated by the vacuum 



energy density, V ~ H 2 Mp. The presence of a no n- vanishing and positive vacuum energy 
density indicates that supergravity is broken and soft masses of order of H are generated [ITT 



The implications of such terms for the charge and color breaking minima and constraints on 
the MSSM parameters have also been recently addressed In minimal supergravity it is 
quite easy to see that such mass terms are generated. In general, the scalar potential in a 



supergravity model is described by a Kahler potential G [13[ and in minimal supergravity 
we define the Kahler potential by 

G = 0*0* + ln|^(0)| 2 (3) 

where <p l represents all scalar fields in the theory and W(<ft) is the superpotential. This 
results in a scalar potential of the form 



V = e 



G 



W i + 6* i W\ 2 l\W\ 2 -3 



(4) 



where W{ = dG/d(f>\ A positive vacuum energy density, V > 0, needed for inflation, breaks 
supergravity and both the exponential and the term enclosed in brackets in flU) must be 
non- vanishing. In fact, the exponential e G is the order parameter for supergravity breaking. 
Included among the 0*'s is the flat direction if, and by inspection of (|]) one finds a mass term 
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e ipip*. A large mass term of this type precludes the possibility of Affleck-Dine baryogenesis. 
It was noted in jlO] that A-D baryogenesis could be revived in non-minimal models in which 
such a mass terms arises with an opposite sign having the effect of driving ip to extremely large 
VEVs. In fact, it was shown in [H|, that in supergravity models which possess a Heisenberg 
symmetry [I5|, including no-scale models of supergravity [15, T7§, supersymmetry breaking 
makes no contribution to scalar masses, leaving supersymmetric flat directions flat at tree- 
level. One-loop corrections in general lift the flat directions, but naturally give small negative 
squared masses ~ —g 2 H 2 /(47r) 2 for all flat directions that do not involve the stop. In the 
context of these theories, we therefore expect initial field values for ip to be quite large. 

In the types of models discussed above, reheating is generally quite inefficient. In models 
in which the inflation is coupled only gravitationally to the observable (gauged) sector, the 
reheating temperature is in fact very low 0, Tr ~ a 2 m? I /Mp ~ 10 5 GeV. However, there has 
been a great deal of recent activity and progress in understanding the quantitative details of 
the process of particle production after cosmological phase transitions and the implications 



of this for reheating after inflation [18, 19, 20, 21, 22], and we would like to briefly discuss the 
relevant aspects of these developments so as to address their implications for our work here. 
One of the key issues is whether or not there is symmetry restoration due to the effects of 
quantum fluctuations It is now clear that the occurrence of this phenomenon is model 
dependent |H% |2"T|, |2"2fl . It has been pointed out [E2J that symmetry restoration effects may 



help localize the fields we are interested in close to the origin. However, these effects can 
occur in the models of interest only if specific couplings between the inflaton, an additional 



scalar field, and the flat-direction field ip are present and assumed to be large WM - In a 
general model, particularly one in which the inflaton has only gravitational interactions with 
other fields, symmetry restoration is not expected to occur. 

In this paper we wish to determine the conditions under which gauge symmetries remain 
broken after inflation due to the presence of additional minima along flat directions. In the 
problem we are interested in, although for sufficiently high reheat temperatures the minimum 
at ip = is a global minimum subsequent to the reheating of the universe after inflation, 
there is another minimum of the potential at ip 7^ which becomes the global minimum 
at lower temperatures and hence later in the evolution of the universe. During the epoch 
when the ip 7^ minimum is the metastable minimum of the potential, it is a relevant and 
appropriate question to ask whether the tunneling rate is such that we will end up at ip = 
due to tunneling. We will examine this rate and hence the possibility of such a tunneling in 
this paper. 



3 



At tree level, the zero temperature potential takes the form 

Vfa Q) = \m 2 (Q)p 2 + -M 4 ">- , (5) 
1 n 

where the mass-squared m 2 (Q) depends on the renormalization group scale Q, and M 
is a new mass scale (such as the Planck mass) which characterizes the strength of non- 
renormalizable interactions due to physics beyond the minimal supersymmetric standard 
model (MSSM). Loop corrections, if included to all orders, would remove the Q dependence 
of V . We will work with the RG-improved tree-level potential, but attempt to optimize the 
choice of Q so as to minimize the contributions of the loop corrections. We treat the case of 
flat directions, where there is no renormalizable interaction term. The lowest value of n which 



is allowed depends on the flat direction of interest P3| ; for the case u r R = s 9 R = b h R = v(Q) 
which we treated earlier, we have n min = 10. For the case u r R = c R = s b R = e# = v(Q), we 
have n m i n = 6. In our numerical investigation of the tunneling rates, we concentrate on the 
case n min = 6, as the rates will be even lower in the case n min = 10. For both flat directions, 
the main contribution to the running of m 2 (Q) comes from SU(3) interactions, which result 
in! 

2 2 2/ ,, Nlrfl f 1 + 3^(M 3 )ln(Q/M 3 )/167r 2 1 

m 2 (Q) = m 2 - —g 2 3 (M 3 )Mi \n(Q/M 3 ) { - / 3 9/ ,, , ^ \ ■ ( 6 ) 



7T^ 



[l + 3(73 2 (M3)ln(g/M3)/87r 2 



where M 3 is the gluino mass at the scale Q = M 3 , and m 2 is the sum of the squark squared 
masses in the flat direction at a scale M 3 . We will ignore the effects of the electroweak gauge 
couplings and the Yukawa couplings, which are sub-dominant. At large enough Q, m 2 (Q) 
becomes negative, and the minimum of the effective potential occurs at ip = v(Q) where 

V (Q) = M (n - 4)/(n - 2) [-m 2 (Q)] 1/(n - 2) . (7) 

The best value of Q to choose (that is, the value which is likely to minimize the higher-loop 
contributions near the minimum) is the solution (if it exists) of Q = v(Q); if no solution 
exists, the apparent minimum is likely to be an invalid artifact of the one- loop approximation. 
Alternatively, we can just choose Q = tp and find the minimum of V(<p, <p). A slightly better 
choice is @ Q = (ip 2 + M 2 ) 1 ^ 2 , which prevents problems at small values of (p, and this is the 
scheme we adopt. 



The temperature-dependent corrections to V(Q, <p) take the form |24 

^(*0 = 2^I> J, dq q 2 \n[l-e ie xp(~[ q 2 + m 2 (p)/T 2 } 1 / 2 
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where T is the temperature, the sum is over all particle species whose tree-level masses 
rrii depend on (p, and 6{ = +1 for a boson and —1 for a fermion. For the flat directions 
we consider, there are a total of 32 particles which couple strongly to the flat direction 
that contribute to the p dependence of Vr(p): 8 gluons with mass g 3 ip, 8 squarks with 
mass (ml + g 2 p 2 ) 1 ^ 2 , 8 gluinos with mass |[(M| + 2g 2 p 2 ) 1 / 2 + M 3 ], and 8 quarks with mass 
i[(M| + 2g 2 (/? 2 ) 1 / 2 — M 3 ]. (Actually, the fermion mass eigenstates are quark/gluino mixtures.) 
The effects of Vr(p) are significant for g 3 p <T, where it can be approximated as 

^)~-|7r 2 T 4 + 2^rV- (9) 

There are also particles which couple to the flat directions weakly and via Yukawa inter- 
actions, and these particles acquire much smaller masses when p is large than those which 
couple strongly to <p. Thus when p > T/g 3 , the number density of the strongly coupled fields 
becomes exponentially suppressed, but the weakly coupled fields can still contribute to Vr, 
and it is therefore these weakly coupled degrees of freedom which are numerically important 
in determining the temperature at which these large scale minima disappear^ These light 
degrees of freedom include bino/quark, hypercharge gauge boson/squark, Higgsino/quark 
and Higgs/squark admixtures, and for sufficiently high temperatures, they contribute an 
amount 

20 i 



to the effective potential, where cq(cio) = y 9 (0) an d {m + ,m_} = {m c , m s }({m u , nib}) 
for n m i n = 6(10). In the n = 10 case, C10 = because there remains an unbroken U(l) 
gauge symmetry which is a linear combination of hypercharge and a color generator. For 
flat directions corresponding to n m i n = 6, the second minimum of Vr (when it exists) occurs 
at v ~ 10 11 GeV, for ~ TeV gluino masses. For n m i n = 10, v ~ 10 15 GeV. 

To calculate the tunneling rate from the (p 7^ minimum to the origin , we need to solve 
the equation 

(5 + ~ r d r ) [V(Q, p) + V T (p)] Q= ^ +Mi)1/2 (11) 

subject to the boundary conditions d r p\ r= o = and p(oo) = v, where v is the location of 
the minimum of V(Q, p) + Vr{p)- The tunneling rate is then proportional to exp(— S3/T), 
where the three-dimensional action is given by 

poo r -, 

S3 = ^J Q drr 2 [l(d r p) 2 + V(Q,p) + V T (p)-V ] J (12) 



-'^We thank Alcssandro Strumia for pointing this out to us. 
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where V$ is a constant which cancels off the potential energy at <p = v. S3 can be computed 
numerically, but to gain intuition we note that at small <p the potential is dominated by 
Vr(p); the desired solution is then 

sinh(2g 3 rr) 

^ r = A 13 

r 

where A is a constant. Near the minimum, the potential can be approximated as 

V(p) ~\m 2 (^-v) 2 (14) 

where m 2 ~ |m 2 (v)| is positive. Near (p — v the desired solution is 

exp(— mr) , , 

r 

where B is another constant. If we now make the drastic approximation that the full potential 
is given by eq.(|9]) for p < p c and by eq.(|14D for <p > p c , with p> c given by the value of <p 
where the two approximate potentials equal each other, then it is possible to solve for A, B, 
and S3 analytically. For T < v and m <C T 2 /v, we find S3/T ~ v 3 /T 3 . Thus we must have 
T near v in order to have a fast enough tunneling rate. 

On the other hand, for T near v , the field-dependent part of Vr dominates the zero- 
temperature potential at p> = v , and for T larger than some critical temperature T crit (mn , M 3 ), 
the second minimum (corresponding to non-zero v ) along the flat direction is removed. In 
fig. 1, we plot contours of T crit (mn, M 3 ), in units of 10 7 GeV, for the n min = 6 flat direction 

— c b r — Sr — &r = V (Q)- We have taken tan/3 = 2. The x-axis specifies the gluino 
mass Ms(M3), while the y-axis labels mq(rhq), where m? = (m 2 — ml R )/3 is the average 
squark mass 2 in the flat direction for rn\ R nig. Above the solid curve, which lies along 
the line rhg ~ O.76M3 for large M3, the zero temperature potential V(ip, Q) (and hence the 
full potential V((p, Q) + Vr) has no second minimum (T cr i t drops off very rapidly to zero 
above the top contour). We see that if the reheating temperature Tr is < 10 7 GeV, the finite 
temperature effective potential after reheating will contain large scale minima for the same 
parameter values as for the zero-temperature effective potential. For larg larger, 
along with the correspondingly larger critical temperatures. In fig. 2, we plot contours of 
T C rit(m n , M 3 ), in units of 10 9 GeV, for the n min = 10 flat direction u r R = s 3 R = b b R = v(Q), 
and we see that for this direction we require larger reheat temperatures, on the order of a 
few xlO 9 GeV, in order to remove the tp 7^ minima. In this case the average squark mass 
rhq is defined by m? = m 2 /3. 

After inflation, the field tp parameterizing the flat direction rolls from p>o, given by (0), 
toward the the nearest minimum of the potential. Because the time-scale for the field 
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to reach this minimum is much shorter than the time for decay of the inflaton and the 
subsequent thermalization of its decay products (assuming, as above, that the inflaton is 
gravitationally coupled to ordinary matter), the evolution of ip is determined by the zero- 
temperature effective potential. The initial conditions at reheating are therefore determined 
by the value of ipo relative to the maximum of V separating the ip = and <p 7^ minima. If 
it is to the left of the barrier, <p will roll to the origin and remain there and we recover the 
scenario described by @. If <po is to the right of the barrier, it will roll to the ip 7^ minimum 
(unless it is so far to the right that it picks up a sufficient velocity so as to carry it over the 
barrier and settle once again at the origin). When non-renormalizable interactions for flat 
directions are included, the maximum value for ipo is limited by the effective mass at large 
ip. For our n = 6 case, the minimum is at roughly ip ~ 10 11 GeV, and the position of the 
barrier depends on m H and M 3 , but is typically much smaller than v, while <^™ ax < 10 14 GeV. 
Thus we expect to be to the right of the barrier, and we have determined that so long as 
ipo < 1.5 x 10 14 GeV, the field will settle in the broken vacuum. For n = 10, the minimum is 
shifted up to about <p ~ 10 15 GeV and is comparable to the maximum value of ipo, while again 
the maximum of the barrier is a much lower scales. Thus if the zero-temperature potential 
for ip has a large-scale minimum, this is where we expect the field to sit at the beginning of 
reheating. If T R < T crit (mn, M 3 ), the position of the minimum of is unaffected by the 
temperature corrections, and after reheating ip will be trapped in the symmetry-breaking 
minimum. 

It is possible that thermal effects can excite ip over the potential barrier and allow nucle- 
ation of bubbles of the symmetric vacuum. In order to determine whether ip remains trapped 
in the ip 7^ minimum, we numerically compute the finite temperature transition rate from 
the ip 7^ minimum to the true minimum (of the finite-temperature potential) at ip — 0. 
For fixed temperature T, we use the full finite-temperature corrections (|8]) to the effective 
potential to compute the bounce solution to (|ll|) and the resulting three-dimensional action 
(|12|). The tunneling probability per volume per unit time is p5 



T ~ T 4 e" 53 / T . (16) 



The fraction of space remaining in the broken phase is / = e p fl26|, where 



and where T is the current temperature of the universe and is the temperature at which 
the two minima become degenerate. For n min = 6 flat directions, is of order 10 6 — 10 7 GeV, 
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while for n m \ n = 10 flat directions Td is of order 10 s — 10 9 GeV; for reheat temperatures less 
than this, there is of course no tunneling. The expression on the right-hand side of (|TT| ) can 
be a gross overestimate, depending on how fast S3 drops off as the temperature falls. As 
P is exponentially sensitive to S3, we can take P = 1 as a critical value, such that P < 1 
implies no nucleation of the symmetric phase. We then find that <p remains trapped in the 
large-scale minimum as long as 

|i >217 + ln fl^ 

For fixed values of (m N , M 3 ), we have looked for values of T R satisfying T R < T CTit (m l} , M 3 ), 
but violating (p!8|). We find that S3 rises very quickly as Tr drops below T crit , and we were 
unable to resolve a single case where (|18D does not hold by better than an order of magnitude. 
That is, Tr must be extremely close to T crit in order for ( |T8| ) to be violated (in which case the 
integral in (|I7|) would need to be evaluated numerically anyway, yielding a much lower value 
for P than given by the right-hand-side of ([T7|)). (In ||, tunneling from tp — to (p — v was 
found to occur over non- negligible regions of the parameter space at small m N , where the 
width of the barrier was small (< 1 TeV); in this case, the barrier widths are much larger). 
We can therefore say that if the reheat temperature Tr is less than the critical temperature 
T cr it (shown in Figure (1)), we expect the field <p parameterizing the flat direction to become 
trapped in the large-scale symmetry-breaking minimum, and to remain there until the large- 
scale minimum becomes the global minimum at low temperatures. This effectively rules out 
regions of (rhq, M3) parameter space where such minima exist, or at the very least implies 
new physics beyond the standard model. 

From fig. 2, we see that typical values of T cr i t for the n m i n = 10 flat direction are about a 
few times 10 9 GeV. We can ask whether reheat temperatures of this magnitude are reason- 
able. There are of course strong limits on the reheat temperature in supersymmetric theories 
due to gravitino production. If the gravitino mass is related to the supersymmetry breaking 
scale (though this need not be the case in no-scale supergravity pTfl ) and m 3 / 2 ~ 1 TeV, then 
the production of gravitinos during the reheat process and their subsequent decay will lead 
to the photodestruction of deuterium and 4 He. This allows one |28| to set a limit on the 
reheat temperature which is comparable to T cr ; t , i.e. Tr < a few xlO 9 GeV, though this limit 
can be evaded for a significantly larger gravitino mass. We note that the gravitino decay 
bound was re-examined in the light of parametric resonance by Allahverdi and Campbell 
2~9fl . They concluded that these bounds cannot be evaded even in the presence of parametric 



resonance. 

In summary, we find that after an inflationary epoch, fields parameterizing flat directions 
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in the SUSY scalar potential will settle into large-scale (> 10 10 GeV) minima, if they exist, 
and will remain there through reheating until the minima become global minima at low 
temperatures. This effectively confirms the qualitative conclusions of albeit with some 
new model-dependence. In particular, we conclude that the SUSY particle mass relationships 
corresponding to the region under the solid line in fig. 1 are inconsistent with inflationary 
models, such as those in no-scale supergravity, in which 1) the inflaton is gravitationally 
coupled to the observable sector, or in any model in which the reheat temperature is less than 
that given by the contours in fig. 1 or fig. 2, depending on which flat direction is populated 
during inflation, and 2) there are no positive 0(H 2 ) contributions to the scalar masses during 
inflation. We note that under these conditions, parametric resonance will not localize the 
fields we are interested in close to the origin. The case of the n m j n = 6 flat directions we have 
considered is the most conservative, in the sense that the regions in the parameter space 
which yield large scale charge and color breaking minima can be made acceptable by lower 
reheat temperatures than in any of the other n m ; n > 6 flat directions involving squarks. The 
^min — 10 flat directions require much higher reheating temperatures. In general, however, 
different flat directions are not necessarily mutually flat, and a large VEV in an n m j n = 6 
direction may preclude large-scale minima in other directions. Conservatism then requires 
that we restrict our conclusions to the n min = 6 case, although there is no reason that 
inflation should populate n min = 6 directions in preference to directions with larger VEVs. 
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Figure Captions 



Fig. 1) Contours of constant T crit , in units of 10 7 GeV, as a function of the gluino mass 
M 3 and the average squark mass rhq = \Jijn^ — ml R )/3, for the n min = 6 flat 
direction u r R = c 9 R = s b R = e R = v(Q). Above the top solid curve, the zero- 
temperature potential (and hence the full potential V(ip, Q) + V T ) has no 
second minimum. 

Fig. 2) Contours of constant T cr i t , in units of 10 9 GeV, as a function of the gluino 
mass M3 and the average squark mass rhq = for the n m j n = 10 flat 

direction u R = s 9 R = b b R = v(Q) . Above the top solid curve, the zero- 
temperature potential (and hence the full potential V(ip,Q) + Vt) has no 
second minimum. 
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